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Non-Abelian anyons — particles whose exchange noncommutatively transforms a system's quantum state — 
are widely sought for the exotic fundamental physics they harbor as well as for quantum computing applications. 
There now exist numerous blueprints for stabilizing the simplest type of non-Abelian anyon, defects binding Ma- 
jorana modes, by judiciously interfacing widely available materials. Following this line of attack, we introduce 
a device fabricated from conventional fractional quantum Hall states and s-wave superconductors that supports 
exotic non-Abelian anyons that bind 'parafermions', which can be viewed as fractionalized Majorana fermions. 
We show that these modes can be experimentally identified (and distinguished from Majoranas) using Joseph- 
son measurements. We also provide a practical recipe for braiding parafermions and show that they give rise to 
non-Abelian statistics. Interestingly, braiding in our setup produces a richer set of topologically protected qubit 
operations when compared to the Majorana case. As a byproduct, we establish a new, experimentally realistic 
Majorana platform in weakly spin-orbit-coupled materials such as GaAs. 

PACS numbers: ()3.67.Lx, 03.65.Vf, 03.67.Pp, ()5.3().Pr 



I. INTRODUCTION 

The search for non-Abelian anyons in condensed matter has 
been a focus of both theoretical and experimental efforts in the 
past decade, driven largely by their potential utility for topo- 
logical quantum computation'. Historically, the first physical 
system thought to host such exotic quasiparticles was a frac- 
tional quantum Hall state at 5/2 filling^. A closely related 
platform for non-Abelian excitations is a chiral p + ip super- 
conductor^. While the existence of intrinsic p + ip supercon- 
ductors in nature remains an open question, Fu and Kane^ pro- 
vided an important insight by showing that heterostructures 
formed by a conventional s-wave superconductor and a topo- 
logical insulator effectively mimic the underlying physics and 
support non-Abelian quasiparticles. Subsequent work showed 
that the topological insulator can be replaced bya spin-orbit- 
coupled semiconductor with Zeeman splitting^. In all of 
these systems, the non-Abelian anyons arise from Majorana 
zero-modes bound to vortices^ Majorana zero-modes are 
also predicted to appear at the ends of ID topological super- 
conducting wires^^ with recent experimental data support- 
ing such a possibility ' ' . Remarkably, in networks of ID wires 
these modes give rise to non-Abelian statistic^^^tSl 

The aforementioned candidate systems feature non-Abelian 
anyons of Ising type, which are not computationally univer- 
sal^. That is, not all unitary quantum gates can be approxi- 
mated by braiding operations alone. While a number ofpro- 
posals to circumvent this problem have been put forwarcC^lE^, 
they are either unrealistic or require non-topological opera- 
tions, hence weakening the main advantage of using anyons — 
topological protection from decoherence. A natural question 
therefore arises as to whether one can find other types of non- 
Abelian anyons in realistic condensed matter systems. While 
there is some hope that the quantum Hall plateau at ~ 12/5 
(and perhaps other fillings) may realize non-Abelian anyons 
with computationally universal braid statistics' '^, experimen- 
tal evidence supporting this idea has yet to appear. 

Following the spirit of recent developments in the pursuit of 



Majorana zero-modes, here we ask whether one can alterna- 
tively "engineer" a composite system supporting more exotic 
non-Abelian anyons. Our search specifically targets quasi-lD 
platforms for such anyons, similar to Majorana wires. At first 
glance, one might be skeptical about the prospects of success 
here given that Fidkowski and Kitaev showed that the only 
non-trivial zero-modes supported by ID electron systems — 
even with strong interactions — are Majorana fermions. We 
will, however, form 'wires' out of edge states of topologi- 
cally nontrivial materials. The 'vacuum' surrounding these 
'wires' is therefore non-trivial, allowing us to circumvent the 
aforementioned restrictions. In particular, we introduce a de- 
vice employing conventional quantum Hall states (such as at 
1/3 filling) and ordinary s-wave superconductors that, remark- 
ably, supports 'parafermionic' zero-modes that in a sense rep- 
resent fractionalized Majorana fermions. These modes gener- 
ate richer non-Abelian braiding statistics that may render them 
better candidates for quantum computation, and can be probed 
via Josephson and tunneling measurements. Our results open 
numerous experimentally relevant directions in the search for 
exotic non-Abelian anyons in condensed matter. 



II. PARAFERMIONS FROM A CLOCK MODEL 

It is instructive to review a toy model discussed recently by 
Fendley^' that supports the parafermion zero-modes that we 
will later realize in a physical electronic system. As a primer, 
let us recall the well-known connection between the transverse 
field Ising model and a spinless p-wave superconductor: the 
former system maps to the latter under a non-local Jordan- 
Wigner transformation that trades the bosonic spin variables 
for fermions. Although the Ising model exhibits only con- 
ventional paramagnetic and ferromagnetic phases, the coiTe- 
sponding superconducting system is far more exotic. Indeed, 
while the paramagnetic phase of the former maps to a triv- 
ial superconducting state in the latter, the ferromagnetic phase 
corresponds to a topological state featuring Majorana zero- 
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FIG. 1 ; Schematic illustration of the parafermion chain Hamiltonian 
in Eq. |4| when (a) J = and (b) /i = 0. In the latter case the ends of 
the chain support 'unpaired' parafermion zero-modes that give rise to 
an TV-fold ground-state degeneracy. 



modes that generate non-Abelian statistics^E^JEil 

Fendley's crucial insight (guided by the identification of 
parafermionic fields in the 2D clock model by Fradkin and 
Kadanof/^ is that one can access still more exotic zero- 
modes by implementing an analogous non-local transforma- 
tion on the generalized A^-state quantum clock model^, 



H = -jJ2ia]<j,+, + H.c.) - hY,{r] + r,). (1) 



Here J > couples neighboring 'spins' ferromagnetic ally, 
/i > is the transverse field, j labels sites of an i-site chain, 
and (Tj , Tj are operators defined on an iV-state Hilbert space 
that satisfy = 1, aj = '^f^'^ similarly for Tj. The 
only non-trivial commutation relation among these operators 



reads OjTj 



r,cr,-e 



2-iTi/N 



When iV = 2 Eq. ([Tji reduces to 



the familiar transverse field Ising model, though the phases 
realized in this special case appear also for general N . For 
example, with J = 0, /i > there exists a unique paramag- 
netic ground state with Tj — +1, while in the J > {), h = Q 
regime an iV-fold degenerate ferromagnetic ground state with 



^ p2mq/N 



emerges (5 = 1,..., N). 



Consider now the non-local transformatiorP3 



a2,=-e'^"'Tja,\{n. (2) 



The properties of (jj , Tj dictate that these new operators satisfy 



af = 1, a] = af -\\nd 



i^sgn{j'-j) 



(3) 



Appendix [A| provides some additional useful relations. For 
N — 2 the operators aj are self-Hermitian, anticommute, and 
square to the identity — hence they are Majorana fermions. At 
larger N, however, they define parafemiion^^^. 
In these variables the Hamiltonian becomes 



"2j 



a2j+i + H.c.) 



3 = i 
N 



h^ie''^ alj_-^a2j + H.c). 



(4) 



The phases of the clock model appear in this representation 
as follows. In the paramagnetic limit with J = the oper- 
ators 'pair up' as sketched in Fig. [TJa) — ai couples to a2, 
couples to 04, and so on down the chain. As shown 
in Appendix [A| one can simultaneously diagonalize this col- 
lection of 'dimers', which take on eigenvalues Q!2j_ia2j = 
_g«Tf ("3-1/2) fpj. jjjtgggj- Here there exists a unique 
ground state with rij = that is fully gapped since exciting 
any of these dimers costs finite energy. More interestingly, 
the ferromagnetic case h = produces the shifted dimeriza- 
tion shown in Fig. [TJb). A bulk gap arises here for the same 
reason, though the ends of the chain now support 'unpaired' 
zero-modes ai and a2L that encode the TV-fold degeneracy of 
the clock model's feiTomagnetic phase {a^j^ai admits TV dis- 
tinct eigenvalues that do not affect the energy). At TV = 2 the 
zero-mode operators ai 2L form the unpaired Majoranas iden- 
tified by Kitaev*^; for TV > 2 they correspond to parafermion 
zero-modeP^that are our main interest here. 



III. PRACTICAL REALIZATION IN QUANTUM WELLS 

The Majorana zero-modes supported by Eq. Q when TV = 
2 are relatively 'easy' to engineer, because in this case the 
operators aj satisfy familiar fermionic anticommutation rela- 
tions [see Eq. ([3]l]. This property allows one to rewrite the 
TV = 2 Hamiltonian in terms of ordinary fermion operators 
Cj = (a2j-i+iQ;2j)/2, yielding a model for a spinlessp-wave 
superconductoi^^ th at can be realized in a variety of experi- 
mental architectures'^^'^^'. Because of the non-standard com- 
mutation relations obeyed by aj with TV > 2, however, devis- 
ing experimental realizations of parafermionic zero-modes is 
substantially more difficult. Our approach is inspired by the 
observation that commutation relations akin to those in Eq. 
Q do occur among physical operators in a familiar system — 
a fractional quantum Hall edge. In particular, for Laughlin 
states at filling factor ly = l/m (m is an odd integer), the 
quasiparticle operators e"^^^'' that create right-moving charge 
e/rn excitations at position x along the edge obey 



(5) 



This suggests that such a system may provide a key building 
block in a device supporting localized parafermion modes. 

A single quantum Hall state is of course insufficient for this 
purpose, since its edge can not be gapped out (and hence one 
can not localize modes of any type at the edge). Thus we con- 
sider the geometry of Fig.|2ja), where two adjacent quantum 
wells, each at filling ly ~ 1 /m, give rise to a pair of coun- 
terpropagating edge states at their interface. These modes 
can acquire a gap via two different mechanisms: first through 
tunneling of electrons across the junction, and second by as- 
sembling electrons from each edge into Cooper pairs. To fa- 
cilitate Cooper pair formation, we will assume that the two 
quantum wells exhibit opposite-sign 5-factors so that the red 
and blue edge states in Fig. [2|a) carry antiparallel spins. (In 
practice one can control the sign of the g-factor by various 
means — see, e.g., Refs. 127128] ) The counterpropagating edge 
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FIG. 2: (a) Experimental architecture realizing parafermion zero- 
modes, (b) Spatial profile for the pairing amplitude /S.{x) and tun- 
neling strength M {x) induced by the superconductors and insulator 
in (a), (c) Schematic dependence of ipix) on the phase difference 
5(f)sc between the superconductors in (a), in the m = 3 case. As 
5(;/>sc winds the larger mismatch between ip{x) on the left and right 
increases the energy until 5(f>sc = 67r. Additional 2ii cycles then 'un- 
twist' ip{x) until the ground state is again accessed at Sipsc ~ 12n. 
Remarkably, this implies that the Josepshon current exhibits 12n pe- 
riodicity in 5(/)sc. 



modes are then similar to those of 2D non-interacting topolog- 
ical insulatorJ^lM] when m = 1 and 'fractional' topological 
insulator^ when m > 1. This allows a pairing gap to open 
at the interface via the proximity effect with ordinary s-wave 
superconductors [green regions in Fig.[2|a)]. A tunneling gap 
meanwhile can originate from spin-orbit-induced backscatter- 
ing between the edge states; the requisite spin-orbit coupling 
can arise either directly from the quantum wells, or from the 
insulator [purple region in Fig. |2|a)] that electrons traverse 
when crossing the interface. 

Of particular interest are the properties of domain walls 
between regions gapped by these different means. One can 



most simply explore this physics by considering the static 
domain structure of Fig. |2|a), in the limit where electrons 
tunnel across the interface only via the central spin-orbit- 
coupled insulator (We assume that spin conservation pre- 
cludes backscattering elsewhere.) We further postulate that 
the potential along the edge is tuned so that low energy right- 
and left-moving e/m quasiparticles carry 'small' momenta. 
In terms of fields (jyj^/L satisfying [(I)r/l{x),(I)r/l{x')] = 
±i^sgn{x — x') and [(I)l{x), (I)r{x')] — i^, such excitations 
are created by operators e*"^^/^ that exhibit commutation re- 
lations of the form in Eq. Q. These properties ensure that 
the electron operators given by 'i/'i?,/L ^ e"""^"/^ obey Fermi 
statistics. Below it will prove useful to write (j)B/L — f ^ 0; 
here p = Oxd/ir is the electron density operator while 



Mx),9{x')] 



i-e{x 



We model the interface with a Hamiltonian H 
wherPl 



Ho 



Hn = 



2ti 



dx [{d.xvf + [d-^ef 



(6) 
(7) 



describes gapless counterpropagating edge modes with 
speed^'' V and Hi encodes couplings induced by the super- 
conductors and spin-orbit-coupled insulator in Fig.|2|a). (For 
simplicity we neglect Coulomb interactions between the edge 
states throughout; also, until specified otherwise we set the 
superconducting phases to zero.) In terms of electron opera- 
tors we have Hi ~ J dx[A{x)^iiil;L + ■M{x)il;]^tpL + H.c], 
where the profiles for the pairing amplitude A(x) and tunnel- 
ing strength M.{x) appear in Fig.|2|b). One can alternatively 
express Hi using ip, 6 variables as 



Hi 



dx[~A{x) cos{2m(p) — A4{x) cos(2m6 



(8) 



Similar models have been studied previously in the m ~ I 
case, where it is well-estab lishe d that each domain wall binds 
a single Majorana fermiorPSSI We will now show that, re- 
markably, for m > 1 the Hamiltonian supports precisely the 
localized parafermion zero-modes that we seek. 

To this end, consider the limit where the induced pairing 
and tunneling terms are sufficiently strong that beneath each 
superconductor (p locks to one of the 2m minima of the first 
term in Eq. ([HI), while under the insulator 9 pins to the minima 
of the second'". Using the coordinates specified in Fig. [2|b) 
one can then write Lp{x < xi) — irnl^^ /m, 6{xi + I < x < 



X2) = nng/m, and ip{x > X2+£) = ttti^ 7m where and 
fig are integer- valued operators. Importantly, Eq. ^ yields 



m 



(9) 



whereas the other integer-valued operators commute. At low 
energies one can focus on the intervals between Xj and Xj + £ 
in Fig. |2|b) where A{x) and Ai{x) simultaneously vanish, 
allowing both ip and 9 to fluctuate. These regions are governed 
by an effective Hamiltonian 



mv \ ^ 

271^^ 

4 — 1 



dx [{dx^f + {djf 



(10) 
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subject to boundary conditions imposed by the adjacent 
gapped regions. As outlined in Appendix |B] the operators 



g"„ vv ' I dx[e 

J Xj 



(11) 



commute with iJgff and thus represent zero-modes bound to 
the domain walls in Fig. |2] Note that aj alters the charge on 
the domain wall by e/m; the first term simply involves e'"^", 
while the other terms similarly add charge e/m. Apart from 
these modes i/gff also supports excitations with a finite-size 
gap of order l/£. We henceforth focus on the ground-state 
sector where such excitations are absent. In this limit the zero- 
mode operators simplify considerably — Appendix [B] demon- 
strates that upon projecting out gapped excitations of iJeff the 



integral in Eq. (Ill collapses to an unimportant constant. One 



then obtains the elegant expressions 



(12) 



that describe the action of the zero-modes within the ground- 
state manifold. 

From Eqs. (|9]l and ( 12 1 it is straightforward to show that 



a^™ = 1, a^aj, = ayaje'^'^<^'-^\ (13) 

Thus the zero-modes bound to our domain walls indeed pro- 
vide a physical realization of the parafermions supported by 
the model discussed in Sec.|ll](with N — 2m). In the present 
context these modes produce a 2m-fold ground state degener- 
acy corresponding to the distinct eigenvalues ofO = a\a2 = 

e*™["»= '""^ '^i/2l.'4ilone can gain intuition for this result by 
noting that since O^hgO — hg + 1, O tunnels between adja- 
cent minima of the cos{2m9) term induced by the spin-orbit- 
coupled insulator The 2m distinct eigenstates of n^^"* — ny^"* 
comprising the ground-state manifold can therefore equiva- 
lently be viewed as linear combinations of the 2m eigenstates 
of fig characterizing the central tunneling-gapped region. 

In the limits studied so far the degeneracy produced by ai 
and a2 is exact. These modes will, however, inevitably hy- 
bridize to some degree — due to quasiparticle tunneling be- 
tween the domain walls — in the more realistic situation where 
ip and 9 are not perfectly pinned by the pairing and tunneling 
terms. For instance, tunneling of e/m quasiparticles between 
the parafermions at lowest order produces a Hamiltonian 



AH = Aa[a2 + H.c. 



'\A\ cos 



(14) 



Due to the gap in the region between the parafermions, the 
coefficient A is exponentially small in Ai and the spacing be- 
tween domain walls. Thus the ground states remain degener- 
ate within exponential accuracy — enabling non-local storage 
of topological qubits in these states. 

Just as for Majorana fermions,'^^ the parafermions in our 
system generate spectacular signatures in Josephson measure- 
ments. To illustrate the physics let us return to the setup of Fig. 



|2]^a) in the limit where </?, are pinned by the superconductors 
and insulator. Suppose that after initializing the system into 
one of the resulting 2m ground states, the tunneling strength 
A4 (x) adiabatically decreases to zero so that the parafermions 
strongly hybridize between the superconductors. [We assume 
that the central region of Fig. |2|a) remains gapped by finite- 
size effects even when Ai{x) ~ 0.] We would like to then 
understand the Josephson current flowing across the junction 
when the superconducting phase on the left side remains zero 
while the phase on the right, S(j)sc, varies. 

The Hamiltonian describing this configuration is H ^ 
Ha + H[, where Hq is again given in Eq. (j?]) while 

H[ — —A j dxcos{2mip) ~ ^ j dxcoa{2m(p — Scpsc)- 

X<Xi X>X2+i 

(15) 

Suppose first that when S(t)sc = we begin in a ground state 
with = ip{x < xi) = and ipy = (p{x > X2 + i) — 
(other initial states are examined below). Upon smoothly in- 
creasing to 2tt our initial state evolves such that ipy = 
5(t>sc/2m to minimize the second term in Eq. ( 15 i. Crucially, 



this cycle raises the system's energy because of the mismatch 
between Lp^ and </?>; the resulting twist of p{x) between 
the superconductors costs energy due to the {dxfY' term in 
Hq. The system returns to its original state only after 6(i>sc 
winds by 47rm, after which one obtains a value = 2ti 
that is physically equivalent to our initial value of 0. Figure 
|2|c) illustrates this physics in the simplest nontrivial case with 
m — i. The Josephson current Io{S(l)sc) oc d{H) /dd(j>sc fol- 
lows from the energy and hence also admits 47rm periodicity. 
(See Appendix [C] for a more quantitative treatment.) 

Interestingly, the current-phase relation depends sensitively 
on the initial state. Consider now the more general situation 
where prior to fusing the parafermions across the junction 
we prepare the system into a ground state characterized by 
= Sn, for some integer dn. One can readily gen- 
eralize the above analysis to show that the current becomes 
Isn{S4>sc) = Ia{S<j>sc + 2TTSn), which indeed differentiates 
all physically distinct values of 6n. Thus the 47rm-periodic 
Josephson effect both provides a definitive signature of the 
parafermions in our setup and enables readout of the quantum 
information they store-. 



IV. PARAFERMION BRAIDING 



The results from Sec. Ill extend straightforwardly to setups 
exhibiting arbitrarily many domain walls separating pairing- 
and tunneling-gapped regions. In particular, 2J\f domain 
walls (numbered sequentially from left to right) localize 
parafermion zero-modes ai.. . 2aA that obey Eq. ( 13 1 and pro- 
duce (2m)-'^ degenerate ground states. Next we show that 
these modes generate non-Abelian statistics. To do so we will 
introduce a practical recipe for transporting domain walls and 
a geometry that permits their meaningful exchange. As we 
will see, compared to the Majorana case braiding enables one 
to perform a richer set of topologically protected operations 
on the ground state manifold. 



To mobilize the domain walls we turn now to the alter- 
nate setup of Fig. [3]^a). Here two quantum wells couple to 
a superconductor and spin-orbit-coupled insulator throughout 
the interface, while gates below control the potential along 
the edges. The Hamiltonian describing the interface is then 
H = Ho 



iJp, with iJg.i given by Eqs. (|7|8|l but now 



with uniform and A. The induced potential ^(x) gener- 
ates the third term, = — J dxii{x)dxd /it (recall that the 
electron density is p = dxd/Ti). We assume A so 

that when ii{x) = a gap arises from inter-edge tunneling. 
Suppose that starting from this regime we adjust the gates 
to raise uniformly. Shifting 9{x) — > 0{x) + -;^x elimi- 
nates the resulting potential term but, crucially, also sends 
M. cos{2m9) — >■ Ai cos{2m6 + 2/ix/w). The spatial oscilla- 
tions render the tunneling term ineffective so that the pairing 
energy A can then dominate. More physically, gating changes 
the momentum carried by low-energy quasiparticles, which 
affects inter-edge tunneling {tp^j^ij^L) but not Cooper pairing 
(V'flV'i)- The gates in Fig. [5]^ a) thereby allow one to dy- 
namically control which regions are tunneling- versus pairing- 
gapped, and hence manipulate domain walls in real time. 

Let us now deform the interface into the 'sack' geome- 
try sketched in Fig. [3|b), where gates are tuned to localize 
parafermions ai . 4 (purple and green respectively denote 
tunneling- and pairing-gapped regions). This setup permits 
exchange of any pair of domain walls shown. We will ana- 
lyze the adiabatic clockwise braid of the domain walls bind- 
ing Qfi and a2 as outlined in Figs. [3|c)-(f); other exchanges 
can be understood similarly. In the figures n\i'' and n^"''' de- 
note, as before, integer operators that describe the pinned </? 
and 9 fields in a given region. Equation (j6|l implies that for 

i > / [n^^Kn"^ \ = iijn/Tr) while for j < j' the oper- 
ators commute. Since we are concerned here only with the 

ground-state sector, it suffices to express ai 



Q!2 



and similarly for a\ 2- 



In the first step [(c)— !-(d)] the right domain wall moves into 
the sack. We model this process with a Hamiltonian 



H. 



-\tj\ cos 
\t\ cos 



H.c.) + [ta'^a'^ + H.c.) 



m \ 



-(2) 



,(3) 



(16) 



The t term above represents charge e/m tunneling between 
the domain walls binding a[ and in Fig. |3jd). Similarly, 
tj reflects e/m tunneling between the inner (red) edge states 
at the junction in Fig.|3jb), with /3 a system-dependent phase. 
Tunneling can also occur between the outer (blue) edge states, 
but only of electrons since the barrier separating the quantum 
wells does not support fractionalized excitations. Electron 
tunneling has no effect on our conclusions, however, and shall 
henceforth be neglected. For details on this important point 
see Appendix[D] The configuration in Fig.[3jc) is described by 
Hc-^d with tj — and i 7^ 0; a'l and a'2 then hybridize while 
ai_2 comprise our initial parafermionic zero-modes. Trans- 
porting the right domain wall modifies the Hamiltonian, re- 
sulting in t J 7^ and t — Q when we arrive at Fig. |3jd). 
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FIG. 3: (a) Setup allowing adiabatic domain wall transport via gat- 
ing, (b) 'Sack' geometry that permits braiding of domain walls be- 
tween tunneling-gapped (purple) and pairing-gapped (green) regions. 
Clockwise exchange of domain walls binding ai and 02 proceeds as 
outlined in (c)-(f). 



Throughout this process ai remains a zero-mode. The second 
zero-mode initially given by 02, however, evolves nontriv- 
ially. In principle one can track the latter by explicitly solving 
Hc^d for general tj, t — a cumbersome task. 

Fortunately, one can deduce the final result far more effi- 
ciently by observing that x = e*2^a2a2 a'l = e*-" '•"4° -* 
commutes with iJc-s-d for any choice of parameters and hence 
is conserved. To see how we first note that at the beginning 
of this step t in Eq. |l6jl pins fSp — n'g^. Consequently 
X projects onto the initial form of the zero-mode of interest, 
X a2- Upon completing this step tj instead imposes the 
energy-minimizing condition 



i(2) 



'(3) 



-(2) 



(17) 



6 



for some /3-dependent integer k. Because x is conserved, 
projecting onto the ground-state manifold yields the properly 



evolved zero-model^ 



a2 «! ■ Notice that this 



zero-mode does not localize to the transported domain wall, 
except in the Majorana case where a'^ — 1. 

Next we transport the domain wall binding ai to arrive at 
Fig.[3|e). This process can be similarly modeled by 



H, 



d— s-e 



{tjala\ + H.c.) + {ia\a2 + H.c.) 



(18) 



where the parameters change from 7^ 0, t = at the 
start of this step to tj — 0, < ^ at the end. In this 
case the operator e^^^^^^^a'2a'i identified above remains 
a zero-mode throughout, whereas now the zero-mode ini- 
tially given by ai evolves. Proceeding as above, we de- 
fine X = e-^^(^-i/2)ai4a; = e*-("l^'+"f '--f '+4^'-^) 
which generically commutes with Hd^e and initially projects 
[using Eq. (fTT])] to ai. When the domain wall moves into the 



sack t pins fig = n'iP ■ The zero-mode formerly described 



(2) I 



k) 



''a'l at 



by Qfi thus evolves to x — 5" e 
the end of this step. 

During the final step of the exchange, where the domain 
walls are transported to the configuration in Fig. [3|f), both 
zero-modes evolve trivially. Let U12 be the operator imple- 
menting the exchange. Comparing Figs.|2|c) and (f), one finds 
that the parafermions transform as 



Ui2aiUl2 

Ul2a2Ul2 



'-(1-k) t 2 



(19) 



Note that a\a2 — the analogue of parity in the Majorana 
case — is preserved here. We would like to now understand 
how the braid transforms the 2m degenerate ground states. 
Assuming ai.2 are the only zero-modes we denote these states 
by \q) where a|a2|(7) = — e^^'^'^^^/^-'lg). Equations ( 19 1 im- 
ply that (up to an overall phase) 



Ui2\q) = e 



(20) 



When additional domain walls are present the braids are 
clearly non-Abelian. As an example, if U23 implements a 
clockwise exchange of the domain walls binding a2 and 

while 



in Fig. gc), then U23U^2a^U\2U, 



-2'Kik/m 



Ui2U2y,aiUl^U\2 = er^^'^/"^a2- Appendix |E| discusses ad- 
ditional braids. In particular, there we show that exchanging 
two pairs of parafermions produces a controlled phase gate 
CP = (C/23C^i2t^34t^23)^ that Can entangle the state of the 
pair ai, a2 with that of the pair 03, 04. Up to an overall phase, 
this operation yields 



CP\q,(^) = e-**f 



(21) 



where q and q' label the eigenvalues of a\a2 and a\a/^ respec- 
tively. Such an entangling gate is unavailable through braid- 
ing of Majorana modes, as can be seen by setting m = 1, 
and indicates a relative increase in computational power for 
parafermions. 



V. DISCUSSION 

In this paper we introduced an experimental setup employ- 
ing conventional quantum Hall edge states to localize ex- 
otic non-Abelian anyons. In the integer quantum Hall case 
(m = 1) our proposal paves the way towards realizing 'Ma- 
jorana wires' in weakly spin-orbit-coupled systems such as 
GaAs. The fractional case (m > 1) constitutes a much more 
important advance, as here our device provides a route to en- 
gineering networks of parafermions that can be moved along 
ID channels. While we focused on a setup with the virtue of 
conceptual clarity, numerous variations on the architecture in- 
troduced here are possible. For example, spin-orbit coupling 
can instead arise from the quantum wells rather than the insu- 
lators in Figs.|2]and|3] The opposite-sign g-factors we invoked 
may also be unnecessary, e.g., if coupling to the parent s-wave 
superconductors does not preserve spin^^ or if one employs 
superconductors with a triplet component^^. We have shown 
that parafermions can be identified — and distinguished from 
Majoranas — via Josephson measurements. Tunneling experi- 
ments provide a less definitive, though perhaps easier, probe 
of parafermions. Indeed, if aj represents a parafermion zero- 
mode then 7j = a™ is a Majorana operator that can couple to 
electrons from a lead; thus parafermions give rise to the same 
quantized zero-bias anomaly as Majoranas^^. Experimental 
verification may also be possible using resonant tunneling of 
e/m charges in a setting where a parafermion mode localizes 
inside a quantum Hall point contact (with tunneling charge 
detectable by, e.g., noise measurement^^. 

In the future it will be interesting to address whether one 
can create even more exotic non-Abelian anyons {e.g., Fi- 
bonacci) using other fractional quantum Hall states beyond 
the Laughlin series considered here. Another worthwhile ex- 
tension of our work would be to explore three-dimensional 
fractional topological insulators with proximity-induced su- 
perconductivity to generalize Fu and Kane's proposal^ for sta- 
bilizing Majoranas. And finally, an important question for 
future applications is whether one can utilize parafermions 
to more readily achieve universal quantum computation. To 
that end, we note that braiding of Majoranas must be supple- 
mented by two additional gates to achieve universality ' . The 
first is a single-qubit phase gate, which together with braid 
transformations completes the set of single-qubit unitary op- 
erations. The second is a multi-qubit entangling gate that may 
be obtained, e.g., by measuring the topological charge of four 
Majoranas. While the first type of gate remains absent in the 
parafermion case, the CP braid described above does gen- 
erate the second type. This gate, combined with arbitrary 
single-'qudit' operations, is universal for quantum computa- 
tion. We leave to future work an investigation of how such 
single-'qudit' operations might be realistically performed. 
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Appendix A: Properties of parafermion operators 

Here we will enumerate several useful properties of 
parafermion operators that follow from the definitions pro- 



vided in SecjllJ Since aj^ — I and = ol^ ^, these oper- 
ators are unitary and exhibit eigenvalues of the form e^'^^il^ 
for integral q. The second equation together with the commu- 
tation relations in Eq. (|3]l imply that 



^sgn(j"'-j) 



(Al) 



Moving ttj' past aj therefore produces the opposite phase fac- 
tor compared to moving aji past aj. Consequently we obtain 
the following commutation relations. 



a.j,ak 
which further imply that 



= {k > i,j or k <i,j), 

(A2) 



alaj,alai 



= 



(A3) 



so long as neither k nor I lie between i and j. 

Equations (A2i and (|A31l demonstrate that as claimed in 



Sec.|ll]one can indeed simultaneously diagonalize each of the 
'dimers' sketched in Figs. flTa) and (b), as well as the combi- 
nation of zero-mode operators a2L'^i case (b). To deduce 
the allowed eigenvalues we note that one can show from the 
properties above that 



(«k)'^ = (-l) 



N+l 



(A4) 



which constrains the eigenvalues of alaj to the form 
—e'lf (9^1/2)^ where q is an integer In the quantum clock 
model context, the eigenvalues of the relevant 'dimer' opera- 
tors can alternatively be found using the relations 



(A5) 




that arise from the non-local transformation specified in Eqs. 

Equations ( |A5| l yield the same eigenvalue spectrum for 
the operators on the left-hand side as noted above since Tj 
and fjj both exhibit non-degenerate eigenvalues e^^^il^ for 
9 = 1,. ..,7V. 



Finally, we consider the case where a\ and a2h represent 
zero-modes and deduce the action of these operators on the 
ground state manifold. Let \q) be a ground state satisfying 
<^\<^2L\'tl — — e*T^'^^^^/^'|g). Using the parafermion com- 
mutation relations one can show that 



(a|a2L)aj = e ' " aj(a|Q!2L) 



(A6) 



for either j = 1 or 2L. This equation implies that q;| 2lI9) 
\q + 1) and ct\.2L\<i) c< |g — 1), where the proportionality 
constants have unit magnitude. One can always fix the relative 
phases of the ground states such that 



a\\q) = \q^\), (l^\q) = \q-\) 
With this convention a2L then acts as follows: 

a\^\q) = -e-'^(''-l/2)|g + l) 



(A7) 



(A8) 



Appendix B: Solution of localized parafermion zero-modes 

This Appendix provides a detailed derivation of the local- 
ized parafermion zero-mode operators quoted in the main text 
[Eq. ([TT])]. Consider again the static domain structure of Fig. 
[2|a). As in the main text we will assume that the Cooper pair- 
ing and inter-edge tunneling terms induced at the interface are 
sufficiently strong that ^p{x) is pinned beneath the supercon- 
ductors while B(x) is pinned beneath the spin-orbit-coupled 
insulator In the black regions of width I in Fig. |2|a), how- 
ever, both fields can fluctuate since the pairing and tunneling 
terms simultaneously vanish there. Our objective is to now 
understand the low-energy properties of these regions, which 
will eventually lead us to the parafermionic zero-mode opera- 
tors of interest. 

For clarity we examine each domain wall separately, begin- 
ning with the one on the left. At low energies this domain wall 
is governed by an effective Hamiltonian 



H = 



27r 



(Bl) 



Minimizing the energy of the adjacent gapped regions re 
quires that these fields also satisfy boundary conditions ip{x - 

xi) = irni^^ /m and 0{x ~ xi + i) = Trhg/m for integer- 
valued operators h^p and hg [recall Eq. (|8|]. Note that h^p 
and hg commute since [ip{xi),9{xi +l)\ = according to 
Eq. ^ 



(1) 



Equation ( B 1 1 can be diagonalized by expanding the ip and 
9 fields as 



ip{x) 
6{x) 



m 
m 



fc=0 



sin Afe(a:) 
^/2k 



cos Afe(a;) 



1 



ak 



(B2) 
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where Aj,(x) 



{2k+l)Tr{x-xi) 
21 



and flfc correspond to conven- 



tional bosonic operators satisfying [afc,a!, 



t 1 - 



[This 



decomposition simuhaneously preserves the commutation re- 
lations amongst (p{x),6{x') and encodes the boundary condi- 



tions specified above.] Inserting Eqs. (B2i into the Hamilto- 
nian yields 



H = ^efc(4afe + l/2) 
(A: + 1/2). 



fc=0 
TTV 



(B3) 
(B4) 



Thus we see that the bosons exhibit a finite-size gap in- 
versely proportional to £. This does not, however, necessarily 
imply that the domain wall admits only gapped excitations. 
We will now demonstrate that, because the operators fi^^^ and 
fig both commute with the Hamiltonian, the domain wall sup- 
ports zero-modes that can be constructed from local operators 
present in our edge theory. 

To do so it is convenient to work with chiral fields (fiR/L 
instead of ip, 9. From Eqs. (B2 1 we have 



± 



2^ [e±*^'=(^)afe + i?.c.] 



E 

fe=0 



V2fc + 1 



(B5) 



which yields the useful relations 

27r 



2-K 



^Lixi+e) = he + q^Rixi+i). (B6) 

m 



One can employ Eqs. ( |B3| l and ( |B5[ ) to show, with the aid of 
various commutator identities and some algebra, that 



(B7) 



Since the right-hand side is a total derivative the above equa- 
tion integrates to 



— ±iv 



(B8) 



Finally, Eqs. (B6 1 and (B8 1 allow one to demonstrate that the 
operator 



xi+i 



+ e~ 



(B9) 



commutes with the Hamiltonian and therefore represents a 
zero-mode of the system. 

Several points are worth emphasizing here. First, ai is 
constructed purely from local e/m quasiparticle operators 
gi-<t>R/L{x) (jj^g represents a physical zero-mode bound to 

the domain wall. To see this explicitly, observe that ni^^ and 

he always appear in ai via e*^"S°' = e^I'^^'^^^+'^^^^i^' and 



e^T^^f = e«[0i?.(^i+f)-0-L(^i+f)]. Second, the expression for 
the zero-mode quoted above is not unique; one can always 
multiply by allowed operators (such as H or e^*^"*) that 
also commute with the Hamiltonian. This freedom — which is 
a generic feature of zero-mode operators — is, however, incon- 
sequential for our purposes. We are concerned here only with 
physics of the ground-state manifold and in operators that cy- 
cle the system across the various ground states, and for this 
purpose the form of ai above suffices. Third, note from Eq. 
( |B5[ ) that the integrand in Eq. ( |B9| l involves only bosonic oper- 
ators flfc, a\ (and not n'^'' or hg). This makes the projection of 

ai into the ground-state manifold with a\ak = very simple; 
upon discarding an overall constant one obtains the result 



ai — e 



(BIO) 



quoted in the main text. 

The right domain wall in Fig. |2|a) can be analyzed very 
similarly. Here the effective low-energy Hamiltonian is given 
by 



H 



mv 

2^r 



X2+i 



dx[{d,^)'' + [djf 



(BID 



where now the fields satisfy boundary conditions 6* (2:2) = 

■nhg/m and Lp[x2 + I) ~ irfvip ' /m. [Since the left and right 
domain walls of Fig. |2|a) are bridged by a single spin-orbit- 
coupled insulator, 9{xi +£) and 6{x2) must be pinned to iden- 
tical values. Thus hg is the same integer-valued operator that 
we introduced above. However, f{xi) and tf{x2 + (-) are 
pinned by different superconductors, which necessitates the 
introduction of distinct integer-valued operators hip' .] Im- 

(2) 

portantly, in this geometry h)p ' and hg no longer commute 

since Eq. (jijl yields [(p{x2 + t)T0(x2)\ ~ {n /m)'^[h^^\hg\ = 
m /m. 

Given our new boundary conditions, the appropriate de- 
composition for Lp and reads 



(fi{x) 



e{x) - 



m 



-4) 



7me 
m 



E 



sin (a:) 



1 



i{ak 



(B12) 



with X'kix) 



as before. 



(2k+l)w(x-x-)) , 1 

i '-^ '- and [ak,a'^,i = uk^k' 

With this expansion our low-energy Hamiltonian once again 
takes the form in Eq. (B3i. One can then follow the steps 

outlined above (taking care to enforce the non-trivial commu- 

(2) 

tation relations between hip ' and hg) to show that the right 
domain wall binds a zero-mode described by an operator 

a2 = e'srK"+"«) dx[e-*^(»$,"+"«)e'^«(-) 



+ e" 



(B13) 



Projecting onto the ground-state manifold yields, up to a con- 
stant. 



a2 



(B14) 
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Appendix C: Parafermion Josephson effect 

In the main text we deduced the qualitative dependence of 
the energy, and hence the Josephson current, on the phase dif- 
ference 6(t)sc between the two superconductors in Fig. |2|a). 
This Appendix explores the physics uncovered there more 
quantitatively. We continue to work in the limit where the 
tunneling strength M [x) vanishes whereas the pairing fields 
in Eq. ( 15 i pin 9 beneath each superconductor The normal 



region between the superconductors can then be described by 
an effective Hamiltonian 




H 



27r 



(CI) 



subject to boundary conditions on f{xi) and 'p{x2 + i) in- 
duced by the neighboring superconductors. 

Because (in contrast to Appendix |B]) the same field is now 
pinned at both endpoints, it is essential that one incorporates 
compactness of ip in what follows; failure to do so yields in- 
correct results for the dependence of the energy on 5(t)sc- We 
will for simplicity set f{xi) = — that is, we fix the eigen- 
value of the operator nl^^ defined earlier to zero without loss 
of generality. At the right boundary, however, we take 



FIG. 4: Energy versus superconducting phase difference 5(j)ac across 
the losephson junction in the m = 3 case. The six curves shown 

(2) 

correspond to the distinct values of n!^, characterizing the pinning 
of ip under the right superconductor, assuming that ip = beneath 
the left superconductor. Provided hip is conserved the energy and 
hence the current are both 12-k periodic in Si^sc- 



with ip{x2 + t) given by Eq. (C2i. The first term in H 



(p{x2 + i) = mod 



(2) 

where hip is the same integer-valued operator introduced pre- 



,(2) 



+ 



2tt 



TT, 27r 



TT, (C2) 



above simply describes gapped excitations with energy ~ 
— ^ — k, which we assume are absent. More interestingly, 
the second term captures the dependence of the energy on the 
superconducting phase difference imposed across the junc- 
tion. 

(2) (2) 

Since [H, hip ] — the eigenvalue of is a conserved 



viously. The right-hand side of Eq. ( |C2[ ) minimizes the pairing 
term in Eq. (15 1 and importantly also restricts (p{x2 + i) to lie 

between -tt and tt for any S(j)sc and n^^^ . Imposing this bound tional' Josephson effects. For a fixed initial value of h^p^ , one 



quantity that can not change under adiabatic evolution of the 
Hamiltonian. It is this crucial property that gives rise to 'frac- 



on the range of ip(x2 + ensures that <^(x) need not exhibit 
any unnecessary twists between x = x\ and X2 + I- 

To diagonaUze the Hamiltonian we decompose (/s, d as fol- 
lows: 



lf(x) = ip{x2 + I) 



X\ 



X2+ i - Xi 



1 



EsinAfc(x) 



(C3) 



sees from Eqs. (C2 1 and (C6 1 that the energy is 47rm periodic 
in S(f>sc, despite the fact that the underlying Hamiltonian — 
recall Eq. ([TSj — clearly exhibits 27r periodicity. [Note that 
here is where compactness of (p is essential. Had we ex- 
pressed ip{x2 + £) in Eq. (C2i without modding by 2tt the 
energy would increase unboundedly with 64>sc, which is obvi- 
ously physically incorrect.] As a concrete illustration, Fig.|4] 
displays the energy £{S(j)sc) versus S(j>sc for the six inequiva- 

^(2 



Oi^) = 0. + :^±'-^ia,+al), (C4) 



k=l 



X2+1 — XI 



and as usual 



are canonical 



where \k{x) 

bosons satisfying [0^,0]^,] = 5k^k'- In Eq. (C4i 6q represents 



the zero-momentum component of 9{x) (note that fc = is 
excluded from both sums above). The boundary conditions 
on (p{x) are clearly obeyed in this representation, while the 
commutation relations between 6 in Eq. ((6]l are also pre- 
served provided (?o and hi^ are conjugate varia bles satisf ying 

[hip ,0q\ — i. Using the decomposition in Eqs. (C3 1 and (C4 1 
one can express the effective Hamiltonian as 



H = ^?fc(4afe + l/2) 

k=l 

mv [^{x2 + i)]'^ 



lent hip ' values in the m — 3 case. 

As mentioned in the main text the Josephson current flow- 
ing across the junction exhibits the same 47rm periodicity as 
the energy. One should, however, bear in mind the following 
caveats that have been raised in the context of the Majorana- 
mediated fractional Josephson effect (see, e.g., Refs. |8l33T l. In 
any experiment the measured current will consist of a 47rm- 
periodic contribution arising from the fused parafermions and 
a conventional 27r-periodic component flowing in parallel. 
(The latter can arise, for example, from the ordinary Joseph- 
son cuiTent that flows directly between the two parent s-wave 
superconductors.) These currents must be disentangled if one 
is to utilize Josephson measurements to read out the qubits en- 
coded by the parafermions. We also note that in practice var- 
ious imperfections — e.g., inelastic processes that change the 
(21 

^{5(t>sc) (C5) value of hip or additional parafermion couplings that spoil 



conservation of h 



(2) 
V ~ 



:an potentially restore 2tt periodicity 



27r X2 + £ — Xi 



(C6) 



of the current. Exploring these subtleties in detail would be 
quite interesting, particularly given the fractionalized nature 
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of the system we are dealing with. 



Appendix E: Controlled-phase gate from sequential braids 



Appendix D: Quasiparticle tunneling at the junction in Fig.[3|b) 



Section IV noted that at the junction in Fig. [3jb) charge 
e/m excitations can tunnel between the right-moving (red) 
edge states, whereas only electrons can tunnel between the 
left-moving (blue) edge states. The distinction between these 
allowed processes is crucial for the outcome of the braid an- 
alyzed there. We now wish to elaborate on this point by ex- 
amining the junction Hamiltonian in greater detail. Let — xq 
and +a:o respectively denote the coordinates of the left and 
right sides of the junction where tunneling takes place. One 
can then model the coupling of right- and left-moving modes 
across the sack by 

H.J = -tiiCos[4)ii{x(^) - (j)R{-XQ) + Pr] 

- tLCOs[m{(j)L{xo) - (j)L{-XQ)) ^ Pl]- (D1) 

Here Ir is the tunneling amplitude for e/m right-moving ex- 
citations, t]^ is the tunneling amplitude for left-moving elec- 
trons, and Pr/l are non-universal phases. Higher-order pro- 
cesses such as pair tunneling can be easily incorporated in 
what follows, but we neglect these for simplicity. 

Our primary interest here is to understand how Hj couples 
parafermion modes. For concreteness let us consider the setup 
in Fig. [sj^d) where zero-modes a2 and a'^ reside on opposite 
sides of the junction. Due to the pinning of (f and 9 near the 
domain walls in the figure, one can replace 



m ^ 



m 



yielding 

Hj = —Ircos 



Im ^ 






-n 








-(2) 



(2; 



(D2) 
(D3) 

) + fiR 
f/3L]. (D4) 



The first line simply corresponds to the {tja\a\ + H.c.) 
parafermion hybridization term in Eq. ( [T6] l. Importantly, the 
second line is a function of the same linear combination of op- 
erators nip -\-nf^ —fi[^^ ""-e'^'' that appears intR above. Con- 
sequently electron tunneling results in a benign coupling of 
the form [t'j(a\a'i)™' + H.c.] that does not change the conclu- 
sions in Sec. |IV] (apart, perhaps, from an unimportant modifi- 
cation of the non-universal parameter k appearing in the braid 
transformation). The factor of tt in the cosine — which reflects 
the fact that describes tunneling of electrons rather than 
charge e/m quasiparticles — underlies this important result. If 
fractionalized quasiparticles could tunnel between right- and 
left-movers, then additional non-local terms would appear in- 
volving parafermions far from the junction. For example, 
tunneling of left-moving e/m quasiparticles would produce a 
term of the form [5t{a2y a'ia2+H .c]. Such couplings, when 
combined with right-moving e/m tunneling, would spoil the 
topological nature of the braids we analyzed but fortunately 
are precluded in our system. 



In this final appendix we consider sequential braids that 
produce a controlled-phase gate. For convenience, we rela- 
bel the ground state basis in terms of 



so that Eq. ( 20 1 yields simply 

C/i2k)'.=e-^^'^|g)',. 



(El) 



(E2) 



Henceforth we drop the label k, assuming that all exchanges 
occur through a single junction characterized by the same k. 

Using this definition, we now examine the effects of more 
complicated braids. As a concrete example consider Fig.|3jb) 
where two pairs of domain walls bind four parafermion zero- 
modes ,4, yielding a ground state space with (2m)'^ states 
(assuming no overall fusion channel constraint). For ease of 
notation let us refer to the domain walls binding aj simply by 
j. One can implement a clockwise exchange of the left and 
right pairs of domain walls via the individual clockwise ex- 
changes of 2 and 3, followed by 3 and 4, 1 and 2, then 2 and 
3 once more, building up the unitary U — J723C^i2C^34C^23- 
In order to determine the effects of this braid, we first fix 
the relative phases of the ground states by defining |p, q)' = 
a\^a]^\Q)' ® |0)', where the first term in the direct product 
denotes the combined state of parafermions ai and a2, while 
the second denotes the combined state of parafermions and 
04. Note that with our conventions we have 



i^{k+m-l/2) 



a\a2\Q)' = — e 



and similarly for 0:3 ^4. 

The total effect of U is to transform 



„ , 27vi(k—l)/m^ 
Ql — > e ^ tta 

27ri(fc-l)/m 



|0)', 



(E3) 



as 



t2 2 



t2 2 
Q!4 — >■ Ckg a20t/^. 



(E4) 



In particular, Ua\a2W = ala4 and Uala^W — a\a2, 
which implies that U\p,q)' — e^'^p? The above rela- 
tions allow one to determine the phases Kp^; upon discarding 
an overall phase we obtain 



u\p,qy 



q,py- 



(E5) 



Double exchange of these two pairs of domain walls thus 
yields, again up to an overall phase. 



u'\p,qy 



(E6) 



corresponding to the CP gate from Eq. (j2T|l upon transform- 
ing back to our original basis. (Note that the junction param- 
eter k cancels out here.) Equation (E61 constitutes a rather 
important result: because the phase factor on the right side 
depends on both p and q, this braid operation can entangle 
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the two registers when acting on a superposition of orthogo- 
nal ground states. This braid distinguishes parafermions from 
Majoranas, since is trivial when m = 1. 

The above braids, together with the {2m)-^ fold degener- 
acy of the ground state manifold, suggest the following set of 
fusion rules for the parafermion modes: 

a0a = Vo ® '^i'l ® • • • ® V'2m-i (E7) 
aigiipj = a (E8) 

V'ji<8V'j2 = V'ji+j2 mod2m- (E9) 



Here xbo is the identity channel and V'i....,2m-i represent the 
distinct quasiparticle types to which pairs of parafermions 
can fuse. (The first line can be intuitively understood from 
our analysis of the parafermionic Josephson effect.) The 
parafermion a has quantum dimension i/2m, while each 
other field has dimension 1. Note that the set of ipj form 
an Abelian sub-algebra consistent with the braid operator U 
found above. For m = 1, these fusion rules reduce to the 
well-known Ising anyon theory. 
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